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ABSTRACT 

Aims. Recent observations of XTE J1739-285 suggest that it contains a neutron star rotating at 1122 Hz. Such rotation 
imposes bounds on the structure of neutron star in XTE J1739-285. These bounds may be used to constrain poorly 
known equation of state of dense matter at densities > 10 15 g cm -3 . 

Methods. One-parameter families of stationary configurations rotating rigidly at 1122 Hz are constructed, using a 
precise 2-D code solving Einstein equations. Hydrostatic equilibrium solutions are tested for stability with respect to 
axi-symmetric perturbations. A set of ten diverse EOSs of neutron stars is considered. Hypothetical strange stars are 
also studied. 

Results. For each EOS, the family of possible neutron star models is limited by the mass shedding limit, corresponding 
to maximum allowed equatorial radius, i? m ax, and by the instability with respect to the axi-symmetric perturbations, 
reached at the minimum allowed equatorial radius, -R m i n - We get R m in — 10 — 13 km, and i? ma x — 16 — 18 km, with 
allowed mass 1.4 — 2.3 Mq. Allowed stars with hyperonic or exotic-phase core are supramassive and have a very narrow 
mass range. Quark star with accreted crust might be allowed, provided such a model is able to reproduce X-ray bursts 
from XTE J1739-285. 

Key words, dense matter - equation of state - stars: neutron - stars: rotation 



o 
o 

(N 

a 

a 
o 

CO 
> 

CN 
O 

Q_i 1. Introduction 



o 

to 

in 



Because of their strong gravity, neutron stars can be very 
rapid rotators; theoretical studies indicate that they could 
rotate at submillisecond periods, i.e., at frequency 



1 /per iod >1000 Hz (e.g., ICook et al.lll994al ISalgado et all 
11994 . During 24 years after its detection, the first mil- 
lisecond pulsar B1 937+21, rotating at / = 641 Hz 
(|Backer et al.l . [l982l) . remained the most rapid one. In 2006, 
a more rapid pulsar J1748-2446ad , rotating at / = 716 Hz, 
was detected (|Hessels et al.l . [200rl ). However, such frequen- 
cies are too low to sign ificantly affect struct ure of neutron 
stars with M > 1M Q (jShapiro et all . 1 19831) . They belong 
to a slow rotation regime, because / is significantly smaller 
than the mass shedding (Keplerian) frequency /k- Under 
the slow rotation regime, effects of rotation on neutron star 
structure are oc (///k) 2 -C 1. To enter the rapid rotation 
regime for M > IMq, one needs submillisecond pulsars 

with / > 1000 Hz 

Very recently iKaaret et all (|2006f ) reported a discov- 
ery of oscillation frequency / = 1122 Hz in an X-ray 
b urst from the X-ray t ransient, XTE J1739-285. According 
to lKaaret et all ([2006) "this oscillation frequency suggests 
that XTE J1739-285 contains the fastest rotating neutron 
star yet found" . The very fact that stable rotation of neu- 
tron star at / > 1000 Hz exists means that instabili- 
ties that could set in at lower / (r-mode instability, other 
gravitational-radiation-reaction instabilities) are effectively 
damped. In the present Letter we derive, using precise 2- 
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D calculations of neutron stars with different EOSs, con- 
straints on neutron star in XTE J1739-285, which result 
from the condition of stable rotation at / = 1122 Hz. 

Neutron star models at / = 1122 Hz are analyzed in 
Sect. [3J where we also derive constraints on neutron star 
parameters resulting from the condition of stable rotation. 
Possibility of strange (quark) star in XTE J1739-285 is 
briefly discussed in Sect. [3l Sect. [4] contains discussion of 
our results and conclusions. 

2. Neutron stars at 1122 Hz 

The stationary configurations of rigidly rotating neu- 
tron stars have been computed in the framework 
of general relativity by solving the Einstein equa- 
tions for stat io nary axi-s ymmetric spacetim e (se e 
iBonazzola et al.l 119931 or iGourgoulhon et all [1999). 
The numerical computations have been performed using 
the Lorene/Codes/Rot_star/rotstar code from the 
LORENE library ( |http : //www . lo rene . obspm . f r ). One- 
parameter families of stationary 2-D configurations were 
calculated for ten EOSs of neutron-star matter, listed 
in Table [1] Stability with respect to the mass-shedding 
from the equator implies that at a given gravitational 
mass M the circumferential equatorial radius R cq should 
be smaller than i? max which corresponds to the mass 
shedding (Keplerian) limit. The value of i? max results 
from the condition that the frequency of a test particle 
at circular equatorial orbit of radius -R max just above the 
equator of the actual rotating star is equal to 1122 Hz. It 
is interesting that the relation between M and R cq at this 
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Fig. 1. Gravitational mass, M, vs. circumferential equa- 
torial radius, R cq , for neutron stars stably rotating at 
/ = 1122 Hz, for ten EOSs labeled as in Tab. [TJ Small- 
radius termination by filled circle: setting-in of instability 
with respect to the axi-symmetric perturbations. Dotted 
segments to the left of the filled circles: configurations un- 
stable with respect to those perturbations. Large-radius ter- 
mination by an open circle: the mass-shedding instability. 
The mass-shedding points are very well fitted by the dashed 
curve R min = 15.52 (M/IAMq) 1 / 3 km. For further expla- 
nation see the text. 



"mass shedding point" is extremely well approximated by 
the formula for the orbital frequency for a test particle 
orbiting at r = R eq in the Schwarzschild space-time created 
by a spherical mass M (which can be replaced by a point 
mass M at r = 0). We denote the orbital frequency of 
such a test particle by f^ w -(M, R cq ). The formula giving 
the locus of points satisfying / s r c b hw (M, R cq ) = 1122 Hz, 
represented by a dash line in Fig. [TJ is 



1 / GM \ 1/2 _ 
27 \R^) ~ ' 



122 Hz 



(1) 



This formula for the Schwarzschild metric coincides with 
that obtained in Newtonian gravity for a point mass M. 
It passes through (or extremely close to) the open circles 
denoting the actual mass shedding (Keplerian) configura- 
tions. This is quite surprising in view of rapid rotation and 
strong flattening of neutron star at the mass-shedding point 
(see Fig. [2]). Equation (UJ implies 



R„ 



15.52 



M 



1.4 M, 



1/3 



km 



(2) 



Out of ten EOSs of neutron star matter, two were chosen 
to represent a soft (BPAL12) and stiff (GN3) extremes of 
the set. These two limiting EOSs should not be considered 
as "realistic" , but they are used just to "bound" the neutron 
star models from the soft and the stiff side. The curves 



Fig. 2. Cross section in the plane passing through the ro- 
tation axis of neutron stars rotating at 1122 Hz, for the 
BGN1H1 EOS ( z > 0) and DH EOS (z < 0). Blank 
area with solid line contour - R e 
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tcq — -ttmaxj hatched area 
-Rmin- The coordinates x and z are defined as 
x = rsin#cos</>, z — rcosO, where r is radial coordinate 
in the space-time metric. 



representing allowed configurations rotating at 1122 Hz are 
plotted in the M — R eq plane in Fig. [JJ 

The EOSs based on most realistic models involving 
only nucleons (FPS, BBB2, DH, APR) result in monotonic 
M(i? cq ) of very characteristic "tilda-like" shape. For these 
EOSs we get i? m in ~ 10 km and i? max w 16 km. The value 
of M(R m i n ) is significantly larger than M(i? max ), the dif- 
ference ranging from 0.5 M@ to 0.8 M Q . The ratio of the 
kinetic to the gravitational energy T/|W| is about 0.09- 
0.11 at the mass shedding limit at -R max and is as low as 
0.03-0.04 at R min . 

Four EOSs are softened at high density either by the 
appearance of hyperons (GNH3, BGN1H1), or a phase 
transition (GMGSm, GMGSp). For them, the range of al- 
lowed masses is very narrow, ~ 0.1 M©, and moreover 
M(R min ) < M(R max ). In spite of (9M/9p c ) /=1122 Hz < 0, 
configurations between i? m i n and i? max are stable with re- 
spect to small axi-symmetric perturbations, because they 
satisfy stability condition 



dM 
dp c 
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(3) 



where the derivative is calculated along a sequence with 
a con stant stellar angular momentum J (|Friedman et all 
1988). Peculiar M(R cq ) dependence is related to the fact 
that configurations belonging to the i? m i n < R oq < Rmax 
are supramassive, i.e., their baryon mass Mb is larger than 
the maximum allowable baryon mass for non-rotating stars, 
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The range of allowed values of R cq is 



not much different from that for purely nucleon EOSs. 
EOSs GMGSp and GMGSm describe nucleon matter with a 
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Table 1. Parameters of neutron stars rotating at / = 1122 
Hz at the i? m i n termination (instability with respect to axi- 
symmetric perturbations) and the i? max termination (mass 
shedding instability) of the M — R cq line, for ten EOSs of 
neutron star matter. 



EOS 


M{Rm in ) 


-Rmin 


^*polc /^"ccj 


T/\W\ 




(Mq) 


(km) 






GN3 a 


2.309 


12.77 


0.813 


0.057 


APR 6 


2.274 


10.47 


0.896 


0.038 


DH C 


2.138 


10.50 


0.887 


0.038 


GNH3 a 


2.119 


12.67 


0.804 


0.055 


BBB2 d 


1.996 


9.928 


0.896 


0.034 


FPS e 


1.872 


9.75 


0.894 


0.033 


BG1H1 / 


1.699 


9.99 


0.876 


0.032 


BPAL12 9 


1.521 


9.64 


0.875 


0.032 


GMHSm' 1 


1.481 


9.65 


0.874 


0.029 


GMGSp' 


1.473 


9.60 


0.875 


0.028 


EOS 


M(iW) 


^max 


f pole/ ^"cq 


T/\W\ 




(Mq) 


(km) 






GN3 a 


2.298 


18.25 


0.564 


0.117 


APR 6 


1.498 


15.93 


0.563 


0.108 


DH C 


1.648 


16.45 


0.560 


0.111 


GNH3 a 


2.163 


17.54 


0.583 


0.103 


BBB2 d 


1.473 


15.85 


0.563 


0.107 


FPS e 


1.393 


15.56 


0.566 


0.104 


BG1H1 / 


1.818 


16.92 


0.568 


0.103 


BPAL12 3 


1.422 


15.61 


0.582 


0.089 


GMHSm' 1 


1.594 


16.19 


0.582 


0.086 


GMGSp' 


1.606 


16.23 


0.581 


0.086 



References fo r the EOSs: ° - iGlendenninel l| 19851 ): b - 
Akmal et all (119981 ): c - iDouchin fc Haensell (|2001l ): d - 
Baldo et al.l (|1997f ). for Paris nucleon-nucleon potential; e - 
Pandharipande fc Ravenhalll (|1989l ) ; f - iBalberg fc Gall (| 19971 ): 
9 - iBombacil (| 19951 ): h - iPons et all (|2000T ). GM+GS model with 
[/ft -" = —130 MeV an d pure normal and kaon condensed phases; 
i - IPons et al.l (|2000l ). GM+GS model with U l t = -130 MeV 
and mixed-phase state between the normal and kaon-condensed 
phases. 



first order phase transition due to kaon condensation. The 
hadronic Lagrangian is the same in both cases. However, 
to get GMGSp we assumed that the phase transition takes 
place between two pure phases and is accompanied by a 
density jump. On the contrary, assuming that the transi- 
tion occurs via a mixed state of two phases, we get EOS 
GMGSm. This last situation prevails when the surface ten- 
sion between the two phases is below a certain critical 
value. The actual value of the surface tension is very poorly 
known, and therefore it is comforting to find that both 
curves are nearly identical. 

3. Quark stars at 1122 Hz 

The possibility that the true ground state of hadronic mat- 
ter is a self-bound pl asma of u, d and s quarks wa s contem- 
plated since 1970s (|Bodmerl [l97ll I Wittenl [l98l . If such 
a state, called strange matter, is the true ground state 
of hadronic matter, qua rk stars built entirely of strange 
quark matter could exist ((Wittcn 1984, Haensel et al.l fl986. 
lAlcock et al| [l986). Could compact star in XTE J1739-285 
be actually a quark star? An accreting quark star would 
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Fig. 3. Gravitational mass, M , vs. circumferential equato- 
rial radius, R eq , for strange stars with crust stably rotating 
at / = 1122 Hz, for several EOS of strange matter. Notation 
as in Fig. Q] EOSs calculated using MIT bag model, with 
a s = 0.2 for all models and m s c 2 = 200, 185, 205 MeV, 
B = 56,45, 67 MeV fm~ 3 for S QM1,2,3, respectively (no- 
tation as in lHaensel et allll986l) . 



have a crust of normal matter and due to instabilities of 
thermonuclear burning of accreted plasma it could be an 
X-ray burster (jHaensel fc Zdunikl Il99ll . iPage fc Cumminj 
2005). The mass of the crust on quark star is at most 
10~ 5 Mq, because the densit y at the crust botto m cannot 
exceed neutron drip density (jAlcock et al.l . ll986l ). 

As we see in Fig. [3j the M(R cq ) curves for quark stars 
terminate at i? max which is rather well approximated by 
Eq. [2] (within 2%). The parameters of non-rotating strange 
stars are subject to the scaling laws with the value of the 
bag constant B (M oc B^ 1 / 2 , R oc B -1 / 2 ) if the strange 
quark mass is scaled as oc .B 1 / 4 (|Zdunikl .r2000). The same 
scaling rules are valid also for strange stars with crust pro- 
vided that the mass of the crust is m uch smaller th an the 
mass of the strange core (see Eq. (7) in lZdunikl (|2002D which 
results in R oc B^ 1 / 2 ). For rotating stars the application 
of these scaling laws requ ires also scaling of the rot ational 
frequency with / oc B 1 / 2 (|Gourgoulhon et aTlll999f ). Thus 
these scaling laws cannot be directly applied to the situ- 
ation presented if Fig. [3] since the models SQMl, SQM2, 
SQM3 rotate with the same frequency 1122 Hz, and the 
condition / 2 = /iy ; B 2 /Bi is not fulfilled. 

In contrast to hadronic EOSs, minimum mass allowed 
for accreting quark star can be very low. However, if we add 
the constraint of > 2.1 ± 0.2 M Q resulti ng from the 

existe nce of a massive slowly rotating pulsar ([Nice et all 
I2005D . we get M(i? max ) > 1.4 M©. Let us notice that the 
values of T/|VT| at i? m ax are much larger than for neutron 
stars, 0.174 and 0.188 for SQMl and SQM2, which makes 
these configuration s secularly unstable with respe ct to the 
triaxial instability (jGondek-Rosinska et aU 120031 ). so that 
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the actual value of i? max could be smaller and M(i? max ) 
could be larger. 

4. Discussion and conclusions 

We constructed neutron star models rotating at 1122 Hz, 
suggested to be the rotation frequency of neutron star in 
XTE J 1739-285. At such frequency, rotation significantly 
affects neutron star structure. The neutron star mass is 
larger then 1.4 Af . The equatorial (circumferential) ra- 
dius is contained between R min = 10 — 13 km and i? max = 
15.52 (AZ/ITA/q) 1 / 3 km = 16 - 18 km. The formula for 
-Rmax(Af) is very precise, despite very strong deformation 
of neutron star at the mass-shedding limit. 

The mass of neutron star in XTE J1739-285 is between 
M(R max ) and M(i? min ). For EOSs assuming purely nu- 
clcon composition of matter, M(R m i n ) is maximum allow- 
able mass at / = 1122 Hz and Af(i? min ) — A/(i? max ) w 
0.5-0.7 M . 

The situation is very different for EOSs with high- 
density softening due to hyperons or due to a phase tran- 
sitions. For such EOSs Af(i? m i n ) is slightly smaller than 
M(ii max ), and the range of allowed masses is very nar- 
row, \M(R m\ n ) - M(R max )\ w 0.1 M G . Moreover, rotat- 
ing neutron stars are then supramassive, with baryon mass 
greater than the maximum allowable baryon mass for non- 
rotating stars, Mb > Af^j^ ax . Their existence requires rapid 
rotation and they collapse into black holes after spinning 
down. However, because of a very narrow range of allowed 
masses, possible presence of such a star in XTE J1739-285 
would be a result of "fine-tuning" in accretion process. This 
could explain the rarity of submillisecond pulsars. However, 
by turning the argument around, this implies that we are 
very lucky to detect an extremely rare system! It should 
be stressed that a measurement o f the mass 2.1 ± 0.2 Af 
of a slowly rotating neutron star (jNice et all . I2005D would 
anyway rule out soft-core neutron stars not only in XTE 
J1739-285, but generally. 

Concerning neutron stars with a phase transition, as 
models of XTE J1739-285, we found that their Af(i? eq ) 
curves are nearly independent on whether phase transition 
occurs between pure phases, or via a mixed-phase state. 

Could compact star in XTE J1739-285 be actually a 
strange star? If so, the strange star should have an ac- 
creted crust of normal matter. It remains to be studied, 
whether the observed X-ray bursts characteristics could 
then be reproduced with such an exotic burster model. In 
contrast to neutron stars, rotation at 1122 Hz alone would 
not require a high AT. However, if we additionally require 
A/max > 2.1±0.2 M Q , then the range of masses is narrower, 
1.4-2.3 M Q . 

Constraints on the dense matter EOS, derived in the 
present Letter, are necessary conditions resulting from the 
existence of neutron star rotating at 1122 Hz, but they are 
not sufficient. Actually, neutron star in XTE J1739-285 ac- 
quired its rapid rotation by accretion. Spin-up by accre- 
tion to 1122 Hz leads to a stronger constraint on the EOS 
than just stable rotation at this freq uency. This problem 
was studied for the 641 Hz pulsar by ICook et al.l (|l994bh . 
Accretion associated with spin can be followed as a line 
in the mass - equatorial radius plane, and 1122 Hz should 
be reached before the star becomes unstable either with 
respect to the mass shedding or the axi-symmetric insta- 
bility. As an example, we studied the accretion-driven spin 



up for the DH EOS. We found that to spin up the star 
by accretion from a Keplerian disk to f — 1122 Hz, one 
needs to accrete 0.37 M for an initial mass Afj = 1.4 Af . 
For Mi = 1.93 Af , accretion of 0.25 Af is needed. For 
higher masses of initially nonrotating star, 1.93 Af < Afj < 
Af max = 2.05 Af , accretion leads to the stellar collapse into 
a black hole before reaching rotation frequency of 1122 Hz. 
Our complete results on the constraints resulting from the 
spin-up by accretion will be presented in a separate paper. 

When this Letter was being completed, the preprint of 
Lavagetto et al. (2006) appeared on the astro-ph server, in 
which neutron star rotating at 1122Hz is used to derive one 
of several constraints on the dense matter EOS. In contrast 
to the present Letter, Lavagetto et al. (2006) use approxi- 
mate "empirical formula" for the mass-shedding frequency 
and get constraints in the mass-radius plane of non- rotating 
neutron stars and strange stars. 
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